Influence of magnetic field offsets on the resistance of magnetic barriers in 

two-dimensional electron gases 
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Magnetic barriers in two-dimensional electron gases are shifted in B space by homogeneous, per- 
pendicular magnetic fields. The magnetoresistance across the barrier shows a characteristic asym- 
metric dip in the regime where the polarity of the homogeneous magnetic field is opposite to that 
one of the magnetic barrier. The measurements are in quantitative agreement with semiclassical 
simulations, which reveal that the magnetoresistance originates from the interplay of snake orbits 
with E x B drift at the edges of the Hall bar and with elastic scattering. 
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I. INTRODUCTION 



The transport properties of two-dimensional electron 
gases (2DEGs) in spatially varying magnetic fields show 
a rich phenomenology, and several interesting effects 
have been reported recently. For example, magnetic 
superlatticesA^ show commensurability oscillations in the 
magnetotransport, while giant magnetoresistance effects 
have been found on one-dimensional magnetic arrays 2 -. 
Also, a variety of single magnetic nanostructures have 
been investigated via their influence on the 2DEG, like 
resistance resonances along magnetic edge states 3 -, the 
experimental realization of magnetic waveguides^, or the 
demonstration of Hall sensing 5 : 6 . A magnetic nanostruc- 
ture of elementary character, known as single magnetic 
barrier - MB, is formed by the highly localized perpen- 
dicular (z-) component Bf z of the magnetic fringe field 
in a 2DEG below the edge of a ferromagnetic film mag- 
netized in transport (x-) directionl&^^^^il 4 *! 5 .. Al- 
ternatively, such a structure can be denned in a 2DEG 
with a graded stepi^. Various aspects of MBs have 
been studied in the past few years. The height of the 
MB can be conveniently tuned by in-plane magnetic 
fields, while their magnetoresistance is interpreted in a 
semiclassical picture ^ 10 ' 11 : 12 ' 13 : 14 ' 15 : 17 , within which elec- 
trons experience a deflection by -B/ z , such that for a 
fixed Fermi energy, the angle of incidence determines 
whether the electrons get transmitted or reflected 7 -^. 
MBs in quantum wires, on the other hand, have at- 
tracted a lot of recent attention from theory. These stud- 
ies have established the potential of MBs as tunable spin 
fflters±Si22£L22^i2S^2L2&2& and predict the presence 
of Fano resonances 3 ^. 

In all experiments reported so far, the resistance of 
MBs has been investigated as a function either of a 
magnetic field applied in x-direction, or of the carrier 
density^^iii^^ii 4 .^^. Therefore, the polarity of the 
total magnetic field in z - direction B z (x) is constant. A 
qualitatively new situation arises when B z (x) changes its 
sign across the MB, since novel types of trajectories with 
a snake orbit character become possible. 
In the present paper, we report the experimental imple- 



mentation of such a magnetic barrier structure. It is 
generated by superimposing a homogeneous perpendic- 
ular magnetic field Bh z to Bf z (x). This way, a MB of 
constant strength and shape is displaced along the B z 
- axis. The most interesting situation arises when B^ z 
partly compensates Bj z (x), such that lines of zero mag- 
netic field along the y-direction exist, which act as guid- 
ing centers for snake orbits. In this range, a pronounced 
dip of asymmetric shape in the magnetoresistance is ob- 
served. We explain this structure by an interplay of the 
conductance enhancement along y direction due to snake 
orbits with both elastic scattering as well as E x B drifts 
at the edges of the Hall bar. Numerical simulations based 
on the Kubo formalism and the Landauer-Biittiker model 
substantiate this interpretation. 

The outline of the paper is as follows: Section II is ded- 
icated to the experimental results. An interpretation of 
the data with the help of simulations is given in Section 
III, followed by a summary in Section IV. 



II. EXPERIMENTS 

The experiments have been performed on a commer- 
cially available GaAs/Al x Gai_ x As - heterostructure^ 1 - 
with a 2DEG 65 nm below the surface. The sample geom- 
etry is shown in Fig. [IJa). The samples were prepared 
by optical lithography, followed by wet chemical etch- 
ing for the definition of the Hall bar, or by metallization 
steps, respectively. Au/Ge Ohmic contacts were defined 
at source and drain contacts and at the voltage probes. 
A dysprosium (Dy) platelet with a thickness of 300 nm 
was deposited on the heterostructure surface by thermal 
evaporation at a base pressure of 2 x 10~ 6 mbar. A Cr/Au 
layer of 150 nm thickness was deposited on top to prevent 
the Dy from oxidizing under ambient conditions. The 
samples were measured in a liquid helium cryostat with 
a variable temperature insert. A rotatable sample stage 
allows orientation of the sample within the x-z plane, see 
Fig. [ijb), with an accuracy of the rotation angle of 0.1 
degrees. The cryostat is equipped with a superconductive 
magnet with a maximum homogeneous magnetic field of 
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FIG. 1: (a) Top view of the schematic sample: a Dy platelet is 
placed on top of a Hall bar with source (S) and drain (D) con- 
tacts as well as voltage probes 1 to 8. (b) Nomenclature used 
for rotated samples, (c) The resistance ARi2{Bhx) of the MB 
(the background resistance of the 2DEG is subtracted), mea- 
sured at a temperature of 2K, and the corresponding magne- 
tization of the Dy film. 



Bh = 8T. A current of I = 100 nA with a frequency 
of 13 Hz was passed from source to drain. Conventional 
four-probe measurements as a function of perpendicu- 
lar magnetic field Bh z were performed, from which an 
electron density of 2.45 x 10 15 m -2 and a mobility of 
22.1m 2 /Vs (corresponding to a Drude scattering time 
of td — 8.4 ps and an elastic mean free path of 1.8 /im) 
were obtained at a temperature of 2 K. The quantum 
scattering time was determined from the envelope of the 
Shubnikov - de Haas oscillations^ 2 - as r q — 1.04 ps to an 
accuracy of 10%. The vanishing of the Hall voltage was 
used to adjust 9 = 0. 

In a previous publication, we have discussed in detail the 
resistance of such a MB as a function of a parallel mag- 
netic field Bhx 1 '' '•■ The fringe field Bf z (x) is given h y 12 i 18 
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see Fig. [2a). Here, zq is the distance of the 2DEG from 
the surface and ho denotes the thickness of the Dy film 
with a magnetization denoted by M(Bh x )- 
The MB was characterized by measuring its resistance 
i?i2 as a function oiBh x at = 0. Furthermore, from the 
fit of Ri2(Bhx) to a semiclassical modet^, the magnetiza- 
tion characteristics of the Dy film at its edge of relevance 
is obtained. Five samples were measured, all with quali- 
tatively identical behavior. Here, we focus on the sample 
with the strongest MB of Bf z ax = B fz (x = 0) = 0.57T 
for Bhx = 8 T, corresponding to a saturation magnetiza- 
tion of n Q M s = 2.1 T. 

With the MB established this way, we superimpose a 
homogeneous, perpendicular magnetic field Bf lz . For this 
purpose, the sample is rotated at B^ = 8 T by small an- 
gles, such that the Dy film is kept at constant magnetiza- 
tion, while a homogeneous magnetic field in z direction 



B hz emerges, see Fig. QJb). In Fig. [SJa), R 18 {B hz ) is 
plotted as a function of Rq^B^z). A linear relation with 
a slope of di?78 /rfi?65 = 1-04 and an offset of 380 SI is 
found. Since both Hall crosses measure the average mag- 
netic field in the probed section, we conclude that the 
MB is not modified by rotations in the range of inter- 
est, namely for \Bh z \ < 1 T, corresponding to |0| < 7.2°. 
Thus, rotating the sample in this range corresponds to a 
displacement of Bf z (x) along the B z - axis, as sketched 
in Fig. [2Kb) . Alternatively, we turned B^ off and set 
= 90° , in which case the MB originates from the rema- 
nence of the Dy film. In that setup, \Bh \ = \Bh z \ < 0.2 T 
can be applied without changing the film magnetization, 
which is kept in-plane by the shape anisotropy. 
In Fig. HKc), the magnetoresistance traces Rvi{Bhz) 
at various temperatures are reproduced. A pronounced 
asymmetric dip of ~ 90S! strength at 2K, located at 
B™ n — — 125 mT, is observed. As the temperature in- 
creases, the dip broadens, both its amplitude and its 
asymmetry are reduced, and the position of the mini- 
mum shifts slightly towards B^ z = 0, but the structure 
remains clearly visible at 32 K. The weak temperature 
dependence of the structure indicates a classical origin. 
In the inset of Fig. [UJc), -B™' 1 is plotted as a func- 
tion of the film magnetization, which is tuned by Bh x - 
B™™ drops approximately linearly as the magnetization 
is reduced. Furthermore, we notice a suppression of the 
Shubnikov - de Haas oscillations for Bh z < 0. 



III. SIMULATION AND INTERPRETATION 

An insightful, yet only qualitative, interpretation 
of the magnetoresistance is possible within the Kubo 
formalism^,, which relates the components <xy, = 
{a;, y}) of the conductivity tensor to the electron velocity 
correlation functions (vi(t)vj(0)) by 



irh 2 



(vi{t)vj(0))dt 



(2) 



Here, the angle brackets denote averaging over all 
trajectories^, and to* = 0.067TO e is the effective elec- 
tron mass in GaAs. Since the formalism yields conduc- 
tivities, its application to a single MB is not meaning- 
ful. Therefore, following Kubrak et ali^, we first con- 
sider a periodic array of MBs, each given by eq. and 
then calculate the resistance per barrier. This technique 
suffers from the difficulty of determining the appropri- 
ate value for the period a, which has been interpreted 
as an active length^ of the MB. For the simulations, 
we used a = 2.5 /im, i.e., such that one period contains 
90% of the integrated magnetic field of a single MB. The 
trajectories of 10 5 electrons with the Fermi energy, in- 
jected in arbitrary directions at y — and arbitrary x 
within one lattice period, are calculated numerically to a 
length of 31 /im. We have assumed small angle scatter- 
ing with Gaussian distributed scattering angles (limited 
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FIG. 2: (a) The Hall resistance i?78 vs. J?65, measured at 2K 
(bold line) and the linear fit to the data (thin line), with a 
slope of 1.04. (b) Sketch of B z (x) = B fz {x) + B hz for B hz =0 
and Bhz = — 125 mT. (c) Measurements of Rw^Bhz)- Inset: 
measured and simulated values of -B™ n as a function of the 
film magnetization. 
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FIG. 3: Components of the conductivity tensor for a periodic 
array of MBs, as simulated within the Kubo model. 



to the range [— tt,tt]) and Poisson distributed times be- 
tween scattering events, with an expectation value equal 
to the measured T q . The full width at half maximum 
of the scattering angle distribution of 0.1577T is given by 
the condition that the simulated value for td has to agree 
with the experimentally determined onei^. 

In Fig. [31 we show the simulated functions <Jij{Bf lz ). 
They have shapes resembling those obtained within the 
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FIG. 4: Main figure: Comparison of the measured MB re- 
sistance as a function of Bhz (full line) with the results of 
simulations based on the Kubo model (triangles) and on the 
Landauer-Biittiker model (circles). The inset shows a sim- 
ulation of R\i(Bhz) for the quasi-ballistic case as obtained 
within the Landauer-Biittiker model. 



Boltzmann model, from which they deviate in three re- 
spects. First, the UijiBhz) traces are displaced in B z - 
direction by —(Bf z ), the average fringe field in one pe- 
riod, and thus centered at vanishing average magnetic 
field, (B z ) = 0. Second, the conductivities are scaled in 
Bhz - and in <7y - direction, with a particularly strong 
enhancement of a yy as compared to a xx . This is a conse- 
quence of electrons moving in snake-type orbits along the 
MB in y-direction, as discussed in detail below. Finally, 
all <7jj are asymmetric around (B z ) = 0. This is most 
clearly visible in the off-diagonal elements, the extremal 
values of which for (B z ) > are significantly smaller than 
those for (B z ) < 0. To compare these simulations with 
the measurements, we use 



. a L 

R\2 = (Pxx ~ Pojr^ + —Po 



(3) 



where L denotes the distance between the voltage probes, 
W the width of the Hall bar and po the resistivity of the 
2DEG in zero magnetic field. p xx is given by 



Px 



' yy 



(4) 



Fig. [4] shows that the simulated magnetoresistance, al- 
beit too low, agrees reasonably with the measured one in 
terms of magnitude and shape of the dip. The asymme- 
try is also visible, although less pronounced than in the 
experiments. The position of the simulated minimum 
corresponds to (B z ) = 0. We note that the simulated 
trace of p yy (Bh z ) remains basically unaffected by the 
MBs, while p xy ((B z )) = —(B z )/ne (not shown). More- 
over, the results for p xx are only weakly dependent on a 
for 2 /zm < a < 3 /im. 

It is well known that periodic, one-dimensional su- 
perstructures (the homogeneous direction being the y- 
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direction) in 2DEGs can increase the diffusion in y- 
direction, which leads to a positive magnetoresistance 
Pxx(B z ) around B z = 02 n 36 i 37 . The situation in our 
experiments is different, since (i) (Bf z ) ^ 0, and (ii) 
Bf z (x) is asymmetric about (B z ) = 0. While the 
first point causes a simple displacement of p X x{Bhz) in 
the Bh z - direction, the second one generates asym- 
metries in p xx (B hz ). In particular, —cr xy a yx (\(B z )\) < 
—cr X y(Jy X {—\{B z )\) which, according to eq. (@}, increases 
P XX (\(B Z )\) as compared to p xx (-\(B z )\). From this 
point of view, our experiments on single MBs are comple- 
mentary to previous studies on periodic one-dimensional 
magnetic superlattices with vanishing average magnetic 
field, where both a positive magnetoresistance around 
B=0 and commensurability oscillations at higher mag- 
netic fields have been found^. Apparently, the positive 
magnetoresistance does not require a periodic structure 
and the related commensurate orbits, but rather emerges 
already from a single unit cell of the superlattice. 

The quantitative differences from the measurements 
have several origins. In the simulated array, trajectories 
that are commensurate with the lattice do exist, but are 
absent in the experiments; the resistances of the MBs 
cannot be added according to Ohm's law as implied by 
eq. ([3]). Also, choosing a value for a remains somewhat 
arbitrary. Furthermore, the simulation does not include 
the edges of the sample. 

These discrepancies can be avoided with the multi- 
terminal Landauer-Buttiker formalis m 38 i 39 . We have 
modeled a sample with sizes as shown in Fig. [I] includ- 
ing source and drain (which was set to zero potential) 
and voltage probes 1-4. The geometry with six contacts 
allows various consistency checks, e.g. simulations of the 
Hall effect outside the MB. The symmetry of the problem 
suggests the equivalence of the set of contacts {S, 1, 2} to 
{D,3,4}. 

We inject 1.8 x 10 5 electrons at each cross-sectional 
line, defined at 2.5 pm inside the contacts S, 1 and 2. 
The electron positions on the injection lines and their 
injection angles are random, while their kinetic energy 
equals the Fermi energy. We solve the equations of mo- 
tion for the electrons and calculate their trajectories until 
they pass the cross-sectional line of any of the six con- 
tacts. The sample edges are simulated as hard walls with 
specular reflectivity. Elastic scattering is included as de- 
scribed above, with scattering times adapted to the ex- 
periments as well as to a quasi-ballistic scenario. Within 
the Landauer-Buttiker formalism, one finds 
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Here the notation 

G J =5Z G ^' i = 1,2,3,4; i = S, 1,2, 3,4, D (6) 

has been used. The individual components of the con- 
ductance tensor are given by 



Gn = —Ni-Tfr j, i = S,l, 2, 3, 4 (7) 

Here, Tji denotes the transmission probability from 
contact i to contact j, and Ni is the number of modes 
in contact i. 

The results of these simulations are reproduced in Fig. 
[4] hi comparison to the experiment. We find good agree- 
ment after a constant resistance of 12 ft is added to the 
simulated trace of R\2{Bh z )- This resistance originates 
from the parabolic magnetoresistivity of the 2DEG as 
a function of Bh x 4 — used to magnetize the Dy film in 
the experiment^ 5 , which at By ix sa 8 T amounts to the 
value taken into account. We emphasize that there is no 
fit parameter in this simulation. The simulated curve re- 
produces all classical aspects of the measurement, namely 
the position of the minimum, the shape of the magnetore- 
sistance dip, and the absolute values of the resistances. 
Similar agreement is obtained for other magnetizations, 
as indicated by the inset in Fig. 2(c), where the sim- 
ulated BT" n is compared with the measured ones. We 
tentatively attribute the remaining deviations to distor- 
tions in the sample not considered in the model, such as 
non-specular scattering and/or finite electric fields at the 
Hall bar edges, inhomogeneities of MB in y direction, the 
spatial variations of Bf x as well as of the electron density 
across the MB, and thermal smearing. 

We proceed by specifying the electron trajectories re- 
sponsible for the enhanced conductance. In Fig. [5l 
four characteristic, ballistic trajectories starting from the 
source contact are shown for a MB with BJ 1 ^ — 0.57 T. 

For Bh z > — 80 mT, the MB is closed, i.e., transmission 
is only possible via E x B drift at the edges of the Hall 
bar. In this regime, three scenarios for electrons entering 
the MB region exist, shown in Figs. [S](a)-(c). Depend- 
ing on the initial condition, i.e. the injection angle of the 
electron into the Hall bar, the electron either moves in 
snake-type orbits along the MB in y-direction, see Figs. 
[5£a,b), or crosses the MB at the lower edge of the Hall 
bar (we denote those orbits as edge orbits, Fig. E£c)). 
At the upper edge, the snake orbits split in two groups. 
They either get reflected or transmitted, depending on 
the point where they hit the upper edge which determines 
the direction of the E x B drift. As a result, a fraction 
of the trajectories that contribute to the enhancement 
of the conductance S yy also participates in the trans- 
mission through the MB, thereby coupling S yy to S xx . 
Transport via edge orbits sets in at B] lz = — 50 mT, and 
its weight increases as Bh z is made more negative. The 
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FIG. 5: Characteristic ballistic trajectories within the dip 
region, (a)-(c), as well as for B z (x) < for all x (d). Bf z (x) 
is indicated in gray scale, while the corresponding B z (x) is 
sketched in the insets. 



snake orbits, on the other hand, are present in the inter- 
val —B™ IX < Bhz < 0, in which B z (x) changes polarity. 
Their weights fluctuate, with an overall tendency to drop 
in the range where edge orbits exist. However, even for 
Bj lz < —Bj 1 ^, the MB has an enhanced conductance in 
y direction, due to orbits of the type shown in Fig. Eld). 
Here, the magnetic field gradient at the MB leads to en- 
hanced backscattering, which may be the reason why the 
Shubnikov-de Haas oscillations are suppressed at nega- 
tive B hz . 

We have calculated R\2{Bh z ) for a quasi-ballistic case 
with a rather high quantum scattering time, namely 
T q = 10.4 ps, which could nonetheless be realized ex- 
perimentally with high mobility heterostructures. The 



result is shown in the inset of Fig. [4] A step rather 
than a dip in i? 12 is obtained at the onset of edge 
orbits across the MB, at which the resistance drops by a 
factor of about 2. Moreover, in the regime where snake 
orbits are present, reproducible resistance fluctuations 
are found, which reflect the varying weight of those 
trajectories. A detailed discussion of these features is 
beyond the scope of this work. The simulations thus 
confirm that the increase of R\2 for Bh z < — 125 mT is 
caused by elastic scattering, an interpretation that is in 
tune with that one obtained within the Kubo model. We 
note that at lower mobilities, scattering between various 
snake orbits smears the simulated resistance fluctuations. 

IV. SUMMARY AND CONCLUSIONS 



The resistance of a magnetic barrier as a function 
of superimposed, homogeneous perpendicular magnetic 
fields has been investigated both experimentally and the- 
oretically. The magnetoresistance shows a characteristic 
asymmetric dip at a homogeneous magnetic field of op- 
posite direction to that one of the barrier. Simulations 
show that electrons can travel along the barrier in snake 
orbits as well as in cyclotron-type orbits with magnetic- 
field gradient induced drifting centers. At the sample 
edges, edge orbits exist due to E x B drifts. They effec- 
tively couple S yy to S xx , thereby reducing the resistance 
of the magnetic barrier. The observed asymmetric resis- 
tance dip is interpreted in terms of the combined effects 
of (i) enhancement of S yy , (ii) transmission via edge or- 
bits and (iii) elastic scattering in the barrier region. Sim- 
ulations based on the multi-terminal Landauer-Buttiker 
model without adjustable parameters result in excellent 
agreement with the experiment. In addition, they predict 
parametric resistance fluctuations in the regime of small 
homogeneous perpendicular magnetic fields, which will 
be interesting to study experimentally in further work. 
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